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SECTION  I 


INTRODUCTION 

Direct  methods  of  obtaining  approximations  to  solutions  of  boundary 
value  problems  often  rely  on  variational  formulations.  Finite  element 
procedures  for  numerical  solution  of  various  engineering  problems  are 
also  based  on  variational  principles.  For  linear  operators  many  inves¬ 
tigators  studied  the  subject  in  an  inner  product  space.  Tonti  [1]  noted 
that  the  self-adjointness  of  an  operator  depends  on  the  given  bilinear 
map.  A  given  operator  may  be  non-self  adjoint  with  respect  to  one  bili¬ 
near  mapping  but  self-adjoint  with  respect  to  another.  Following  Tonti, 
Magri  [2]  showed  that  for  every  linear  operator  there  is  an  infinity  of 
bilinear  mappings  with  respect  to  which  it  is  self-adjoint.  Guntin 
[3,4]  used  convolution  product  as  a  bilinear  mapping  for  initial -bound¬ 
ary  value  problems.  Sandhu  and  Pister  [5,6]  extended  the  application  of 
this  bilinear  mapping  to  linear  coupled  initial-boundary  value  problems 
using  a  generalization  of  Mikhlin's  basic  theorem  [7]. 

Mikhlin  [7]  assumed  homogeneous  boundary  conditions  in  stating  the 
variational  principle  so  that  a  self-adjoint  operator  was  symmetric 
[10J.  The  conventional  procedure  to  treat  nonhomogeneous  boundary  condi¬ 
tions  has  been  to  find  the  particular  solution  for  the  nonhomogeneous 
boundary  conditions  and  use  change  of  the  field  variable  to  homogenize 
the  boundary  conditions.  This  approach,  though  theoretically  elegant,  is 


cumbersome  in  implementation  of  direct  methods.  Gurtin  [3,4]  introduced 
nonhomogeneous  boundary  terms  explicitly  into  the  governing  function. 
Nickell  and  Sackman  [8]  and  Sandhu  and  Pister  [6]  followed  Gurtin's  ap¬ 
proach.  In  the  context  of  application  of  finite  element  method,  Prager 
[9]  included,  in  the  variational  formulation,  jump  discontinuities  which 
may  exist  across  interelement  boundaries.  This  development  assumes  spe¬ 
cial  significance  in  case  continuity  of  approximating  functions  cannot 
be  ensured  up  to  the  desired  degree.  Sandhu  and  Salaam  L 10 J  examined 
the  general  case  of  linear  operators  with  nonhomogeneous  boundary  condi¬ 
tions  and  internal  jump  discontinuities  based  on  extension  of  Mikhlin's 
theorem.  By  introducing  the  concept  of  boundary  operators  consistent 
with  the  field  operators,  a  systematic  procedure  to  obtain  variational 
principles  for  linear  coupled  problem  was  developed.  In  line  with  this 
approach,  Sandhu  [11]  presented  a  comprehensive  treatment  on  the  varia¬ 
tional  principles  for  soil  consolidation  using  convolution  product  as 
the  bilinear  mapping. 

Ghaboussi  and  Wilson  [12]  derived  a  variational  principle  for  the  dy¬ 
namic  analysis  of  saturated  porous  elastic  soil.  Biot's  equations  of  mo¬ 
tion  [14]  were  restated  in  integral  form  through  Laplace  transformation 
followed  by  rearrangement  of  terms  and  inversion.  The  general  procedure 
followed  Sandhu  and  Pister  [13,14].  Herrera  and  Bielak  L 1 5 J  pointed  out 
that  the  equivalent  variational  formulation  could  be  written  without 
transforming  the  field  equations  and,  instead,  using  Tonti’s  approach 
[16J.  Ghaboussi's  [12]  treatment  of  the  boundary  conditions  was  inconi- 


the  boundary  conditions  in  a  consistent  fashion  and  to  develop  a  system¬ 
atic  procedure  for  variational  formulation  governing  dynamics  of  fluid- 
saturated  linear  porous  media.  Both  direct  and  complementary  formulation 
of  field  equation  are  studied  following  Sandhu's  approach.  In  chapter 
II,  some  mathematical  concepts  and  definitions  basic  to  the  development 
of  variational  principles  are  introduced.  Biot's  field  equations  for  dy¬ 
namics  of  fluid-saturated  linear  elastic  porous  media  are  given  in  Chap¬ 
ter  III.  Integral  form  of  the  equations  is  the  same  as  used  by  Ghaboussi 
and  Wilson  [12],  In  Chapter  IV,  general  variational  principles  for  the 
problem  are  developed  following  Sandhu  and  Salaam  [10J  and  Sandhu  [11]. 
The  governing  functional  for  the  operator  equations  explicitly  includes 
the  initial  conditions,  the  nonhomogeneous  boundary  conditions  as  well 
as  any  internal  jump  discontinuities.  As  an  alternative  procedure,  com¬ 
plementary  formulation  of  the  problem  is  presented.  Extended  variation¬ 
al  formulations  based  on  self-adjointness  of  the  operator  matrix  are  in¬ 
troduced  along  with  several  specializations.  These  should  be  useful  as 
the  starting  points  for  alternative  approachs  to  finite  element  formula- 


SECTION  II 


MATHEMATICAL  PRELIMINARIES 


2.1  BOUNDARY  VALUE  PROBLEM 

Consider  a  boundary  value  problem 

A  ( u )  ~  f  on  R 

(1) 

C(u)  =  g  on  SR 

(2) 

where  R  Is  an  open  connected  region  in  an  Euclidean  space.  bR  is  the 
boundary  of  R  and  R  its  closure.  We  suppose  that  the  field  operator  A 
and  the  boundary  operator  C  are  bounded  and  defined  such  that 

A  :  wR  - >  VR 

C  :  WbR  >  VbR 

VR,  V^r  are  linear  vector  spaces  defined  on  the  regions  indicated  by  the 
subscripts  and  WR,  W^R  are  dense  subsets  in  VR,  V^r,  respectively.  Thro¬ 
ughout,  A  and  C  are  assumed  to  be  linear  so  that 

A  (a  u  +/3v)  =  a  A ( u )  +  3 A ( v )  V  u,v£WR  (3) 

and  C  (a  u  +/9v)  =  a  C(u)  +/9C(v)  V  u,v£W&R 

-  4  - 


(4) 


where  a,  /3  are  arbitrary  scalars. 


2.2  BILINEAR  MAPPING 

A  bilinear  mapping  B  :  WxV  - >  S,  where  W,  V,  S  are  linear  vector 

spaces,  for  given  w£W,  v£V,  is  defined  as  a  function  to  assign  an  ele¬ 
ment  in  S  corresponding  to  an  ordered  pair  (w,v).  B  is  said  to  be  bili- 
nea  r  i f 

B(a  Wj  +/3w2,  v)  =  a  B(w.  ,v)  +  /3B(w2,v)  (5) 

B ( w,  a  v1  +$v2)  =  aB(w,Vj)  +  £B(w,v2)  (6) 

where  a,  /3  are  scalars.  We  shall  use  the  notation 

Br(w,v)  =  <w,v>R  (7) 

B  is  said  to  be  nondegenerate  if 

<w,v>R  =  0  V  w€W  if  and  only  if  v=  0  (8) 

For  W  =  V,  B  is  symmetric  if 


<w,v>R  =  <v,w>R 


(9) 


■r^y' 


2.3  SELF-ADJOINT  OPERATOR 


Let  A:  V  — >  W  be  an  operator  on  the  linear  vector  space  V  defined  on 

★ 

spatial  region  R.  Operator  A  is  said  to  be  adjoint  of  A  with  respect  to 
a  bilinear  mapping  <  ,  >R  :  W  x  W  — >  S  if 

<w,Av>r  =  <v,A*w>r  +  DaR(v,w)  (10) 

for  all  weW  and  v£V.  Here,  D^R(v,w)  represents  quantities  associated 

★ 

with  the  boundary  bR  of  R.  If  A  =  A  ,  then  A  is  said  to  be  self-adjoint. 
In  particular,  a  self-adjoint  operator  A  on  V  is  symmetric  with  respect 
to  the  bilinear  mapping  if  V  =  W  and 

<w,Av>r  =  <v,Aw>r  (11) 

2.4  GATEAUX  DIFFERENTIAL  OF  A  FUNCTION 

The  Gateaux  differential  of  a  continuous  function  F  :  V  — >  S  is  de¬ 
fined  as 

1 

A  F(u)  =  lim  —  [F(u+  Av)  -  F(u)]  (12) 

v  A->0  A 

provided  the  limit  exists,  v  is  referred  to  as  the  'path'  and  A  is  a 
scalar.  We  note  that  for  u,v£V,  u+Av£V.  Equation  (12)  can  be  equiva¬ 
lently  written  as 

A..F(u)  =  —  F(u+  A  v  )| 


(13) 


2.5  BASIC  VARIATIONAL  PRINCIPLE 


For  the  boundary  value  problem  given  by  (1)  with  homogeneous  boundary 
condition,  Mikhlin  [7]  showed  that  for  self-adjoint,  positive  definite 
operator  A,  the  unique  solution  uQ  minimizes  the  functional 

fl(u)  =  <Au,u>R  -  2<u,f>R  (14) 

where  <  ,  >R  denotes  inner  product  over  the  separable  space  of  square 
integrable  functions.  Conversely,  uQ  which  minimizes  the  functional  (14) 
is  the  solution  of  the  problem  (1). 

Taking  Gateaux  differential  of  (14), 

1 

AQ(u)  =  lim  — [<A(u  +  Av),u  +  Av>  -  2<u  +  Av,f>  -  <Au,u>  +  2<u,f>] 
v  A-0  A 

=  <Au,v>  +  <Av,u>  -2<v,f> 

=  2<v,Au-f>  =  0  (15) 

In  writing  (15)  we  have  only  assumed  linearity  and  self-adjointness  of  A 
with  respect  to  the  bilinear  mapping  and  the  symmetry  of  the  bilinear 
mapping.  The  Gateaux  differential  evidently  vanishes  at  the  solution  uQ 
such  that  AuQ  -  f  =  0.  For  the  vanishing  of  the  Gateaux  differential  at 
u  =  uQ  to  imply  AuQ-f  =  0,  the  bilinear  mapping  has  to  be  nondegenerate. 
To  prove  the  minimization  property,  the  bilinear  mapping  has  to  be  into 
the  real  line  and  the  operator  must  be  positive.  However,  in  general,  it 
is  only  necessary  to  use  vanishing  of  the  Gateaux  differential  as  equiv¬ 
alent  to  (1)  being  satisfied. 
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SECTION  III 

FIELD  EQUATIONS  FOR  DYNAMICS  OF  FLUID-SATURATED  POROUS 

SOILS 


I 

9 


$ 

s.' 


3.1  PRELIMINARIES 

In  this  section,  field  equations  for  dynamics  of  fluid-saturated  lin¬ 
ear  elastic  porous  media  are  stated  following  Biot  [13,14].  The  domain 
of  definition  of  all  functions  is  the  cartesian  product  Rx[0,»),  where 
R  is  the  closure  of  the  spatial  region  R  and  [0,«)  is  the  positive  time 
interval.  The  soil  skeleton  is  assumed  to  be  linear  elastic.  Integral 
form  of  the  field  equations  is  obtained  by  Laplace  transformation  of  the 
Biot's  equatins  of  motion  followed  by  inversion.  Throughout,  standard 
indicial  notation  is  used.  The  Latin  indices  take  on  range  of  values 
1,2,3  and  summation  on  repeated  indices  is  implied.  A  superposed  dot  in¬ 
dicates  time  derivative. 

3.2  DIFFERENTIAL  FORM  OF  FIELD  EQUATIONS 
(a)  Dynamic  Equilibrium 

Biot's  equations  of  motion  for  the  binary  mixture  of  fluid  and  solid 
are  L14J; 

■Clj.J  tPt>i  ■  p51  +  p2Sf  1,61 


7T  ,  +  Pob-  =  P.j'u.  +  (Po/f)w,  +  (l/k)w. 


(17) 


•vi'ivr 


where  'g.y  b^ ,  u^ ,  w^  are,  respectively,  components  of  the  total  stress 
tensor,  the  body  force  vector  per  unit  volume  of  the  mixture,  the  solid 


displacement  vector  and  the  relative  displacement  of  the  fluid  with  res¬ 


pect  to  the  solid.  TT  is  the  pore  pressure.  P  is  the  mixture  density  and 


P 2  is  the  mass  of  fluid  per  unit  volume  of  the  mixture. 


(b)  Kinematics 


For  small  displacement,  the  strain-displacement  relations  are; 


e1j  =  u(i,j) 


*1/2 


£  =  w .  * 

5  i,1 


where  e. .  are  components  of  the  symmetric  strain  tensor  of  solid  and  £ 

•  J 

is  the  rate  of  volume  change  of  the  fluid  per  unit  volume  of  the  solid. 


(c)  Constitutive  Equations 


For  linear  elastic  fluid-saturated  soil,  Biot  [14]  proposed  the  con¬ 


stitutive  equations 


=  Eijkl  ekl  *  aM  5ij  ,a5kl  ekl 


V  -  M  ,a51j  efj  +f> 


The  inverse  relationships  are; 


ei  j  =  Ci  jkl  (  5kl  "  a7r5kl ) 


f  =  7r(l/M  +  a  C.jkl  5k1  5..)  -  aC.jkl  5R1  £ 


K 


R 


K* 


1 


.4 


v -ic.  v  vv-* . C\C\-C 


Here,  and  C^^  are,  respectively,  components  of  the  elasticity 
and  compliance  tensors  of  the  elastic  solid,  a  is  the  compressibility 
of  the  solid  and  M  that  of  the  fluid. 


{ d )  Boundary  Conditions 


u4(x,t)  =  u,(x,t)  on  S  x  [0, x ) 

1  -  1  -  i 

w.(x,t)  =  w.(x,t)  on  S  x  [0, x ) 

1  ~  1  2 

X.-fx.tln.  =  T  •  (x,t)  *  T  -  ( x ,  t )  on  S  x  [0,ac) 

’  J  J  •  •  O 


^(x.tjn.  =  7T.(x,t)  on  S  x  [0,x) 
1  1  '  4 


The  first  two  in  (24)  are  the  displacement  boundary  conditions  and  the 
last  two  the  traction  boundary  conditions,  n^  are  the  components  of  the 

A 

unit  outward  normal  to  OR  and  T..  are  components  of  the  prescribed  tract 
ion  in  i-th  direction.  Each  of  {$1,83}  and  |S2,S4]  consists  of  disjoint 
complementary  subsets  of  OR. 


(e)  Ini tial  Conditions 

The  initial  conditions  for  the  problem  are 


u  (x,0 )  =  u  (x) 


u  (x,0)  =  u  ( x ) 


w  (x,0)  =  w  (x) 


SO 


m 


x 


The  equations  (16)  through  (25)  completely  define  the  initial-boundary 
value  problem  of  small  motion  of  fluid-saturated  porous  media. 


3.3  INTEGRAL  FORM  OF  THE  FIELD  EQUATIONS 

For  development  of  variational  principles,  we  need  to  rewrite  the 
field  equations  in  the  form  of  convolution  product  so  that  the  time  de¬ 
rivatives  are  avoided.  This  can  be  done  through  applying  Laplace  trans¬ 
form  and  taking  inverse  after  appropriate  rearrangement.  Following  the 
procedure  originally  suggested  by  Gurtin  [3,4],  Ghaboussi  and  Wilson 
[12]  presented  the  following  results  for  the  Equations  (16)  through 
(25). 


(a)  Dynamic  Equilibrium 

Laplace  transformation  of  (16)  and  (17)  followed  by  inversion  gives; 


**  rij.j  *  Fi  -  -  P2“1  ■  0 

t*  n,  -  l  *  ( 1  /k )W|  +  G,  -  P2  u,  -  (P2/f)  wf 


where 


F^  =  t*  Pb-  +  P[tu1-(0)  +  u.j(0)]  +  p^  [tw.j(0)  +  w^O)]  (28) 

Gi  =  l*  P2bi  +  P2  +  ui(0)]  +  (P2/f)  Ctw^O)  +  wl(0)] 


+  (1/k)  tw. (0 ) 


(29) 
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Here,  symbol  *  denotes  the  convolution  product  defined  as 


f  *  g  -  f  f(  r  )  g(t-  r  )  dT 
•'o 


(30) 


which  satisfies  the  commutative  law,  the  associative  law  and  the  dis¬ 
tributive  law. 


(b)  Kinematics 

Equations  (18)  and  (19)  need  to  be  restated  in  the  form  [11,12] 


t*  e..  =  (1/2)  t*  (u,  .  +  u,  ,) 

I  J  '  9  J  J  9  I 


t*  £  =  t*  w 


1,1 


(31) 

(32) 


{ c )  Constitutive  Equations 

Equations  (20)  to  (23)  must  be  restated  so  that  the  constitutive  re¬ 
lations  show  the  dependence  of  quantities  appearing  in  the  equilibrium 
equations  upon  corresponding  kinematic  quantities  in  them. 


t*  rjj.  =  t*  Eijkl  ekl  +  t*  a  M  5ij  (a.Sklekl  +  £)  (33) 

t*7r  =  t*  M  (o(5. .  e.  •  +  £)  (34) 

I  J  '  J 

**  etj  ■  ‘*c1Jkl  <tki  '  «**ki>  (35) 

t*  f  -  *  a2  CjJk,  5k,  S(j)  -  t*  ac1jkl  «kl  t,j  W 


i 


SECTION  IV 


VARIATIONAL  PRINCIPLES 

4.1  PRELIMINARIES 

To  transform  the  coupled  initial-boundary  .value  problem  of  wave  equa¬ 
tions  for  the  fluid-saturated  porous  media  into  an  equivalent  variation¬ 
al  problem,  the  set  of  field  variables  are  regarded  as  an  n-tuple  in  the 
admissible  space  whose  elements  are  defined  in  Rx[0,°°).  A  solution  of 
the  mixed  problem  is,  then,  an  admissible  state  of  the  field  variables 
which  satisfies  the  field  equation,  the  initial  conditions  and  the 
boundary  conditions  to  the  problem.  The  linear  vector  space  W  consisting 
of  all  admissible  states  is  referred  to  as  the  product  space,  i.e. 

W  =  W,  x  W9  x  . .  x  W„  (37) 

i  c  n 

where  w.  is  an  subspace  whose  elements  represent  the  admissible  state 
for  an  specific  field  variable  u... 

Consider  operator  equations  of  the  coupled  boundary  value  problem 


A  ( u )  -  f  on  R 

(38) 

C(u)  =  g  on  OR 

(39) 

-  13  - 


I 


I 

1 


in  which  u  is  the  vector  of  the  field  variables  and  A  is  an  operator 
matrix  which  is  self-adjoint  in  the  sense  of  an  appropriate  bilinear 
mapping.  Each  element  of  the  operator  matrix  A  may  be  viewed  as  a  trans¬ 
formation 


P. 


(40) 


I 


where  M^,  are,  respectively,  the  domain  and  the  range  of  A^j  which 
are  both  the  linear  vector  spaces  on  R.  The  explicit  form  of  the  field 
equations  is 


of  Mikhlin's  theorem,  the  governing  function  for  the  operator  equation 
(38)  and  (39)  is  defined  as 


J 

For  the  present  problem, 
tin  viz. 


-  2f(>R  ♦  <u,,CjjUj  -  2s, >w)  («) 

we  use  the  bilinear  mapping  introduced  by  Gur- 


(46) 


4.2  CONSISTENT  BOUNDARY  CONDITIONS  AND  INTERNAL  DISCONTINUITY 

Sandhu  [10]  pointed  out  that  appropriate  boundary  terms  should  be  in¬ 
cluded  in  the  governing  function  even  if  they  are  homogeneous.  This  is 
important  for  certain  approximation  procedures,  e.g.  the  finite  element 
method,  where  the  functions  of  limited  smoothness  are  used.  The  bound¬ 
ary  operators  must  be  in  a  form  consistent  [10]  with  the  field  operator. 
Consider  the  boundary  value  problem  of  multi-variables  given  by  (41)  and 
(42).  Referring  to  (44),  Sandhu  [17]  defined  consistency  of  boundary  op¬ 
erators  with  the  field  operators  to  be  the  property; 


n  n 

D5R(ui,uj)  =  <vi’  SCijuj>0R  '  2/Uj,Cji VbR  *  1=1>2’“,,n  (47 

j  j 


In  seeking  approximation  to  the  exact  solution  by  the  finite  element 
method,  the  function  space  with  limited  smoothness  over  the  entire  do- 


main  is  sometimes  used.  In  order  to  properly  handle  this  limited 
smoothness  problem  in  the  variational  formulation,  Sandhu  [11]  intro¬ 
duced  internal  discontinuity  conditions  in  the  form; 

(Cu) 1  =  g  on  5Ri  (48) 

where  a  prime  denotes  the  internal  jump  discontinuity  along  element  bou- 
dary  oR^  embedded  in  the  region  R.  Sandhu  and  Salaam  [10]  and  Sandhu  [11] 
showed  that  this  condition  can  be  included  explicitly  in  the  governing 
function. 

4.3  VARIATIONAL  PRINCIPLES  FOR  DYNAMICS  OF  FLUID-SATURATED  SOILS 
4.3.1  Field  Equations 

Equations  (26)  through  (36)  in  self-adjoint  matrix  form  are; 

A(u)  =  f  on  R  x  [0, w  )  (49) 

Here, 

p  P2  0  -L  0  o' 

P,  Po/f+l*(l/k )  -t*-  0  0  0 

*■  1  am 

0  t*—  0  0  0  -t* 

dm 

A  =  (50) 

L  0  0  0  -t*  0 

0  0  0  -t*  P  t*cM6i  j 

0  0  -t*  0  t*  a  M  fik]  t*M 


I 


Consistent  form  of  the  internal  jump  discontinuities  is 


-t* 

-t* 


t* 


t* 


<u,nj)' 

=  -t*  (g  ) • n • 
1  1  J 

on  S^.  x  [0,oo  ) 

(wini>' 

=  -t*  g 

2 

on  $2.j  x  [0,x  ) 

(7^.)' 

=  t*  g  n. 

3  1 

on  S3l-  x  [0,»  ) 

) 1  =  t*  g  n . 

4  1 

8 

o 

1 _ 1 

X 

oo 

c 

o 

(55) 


Here,  surfaces  S^.,  $2^,  S3l-  and  S4l-  are  embedded  in  the  interior  of  R. 
Operators  in  the  self-adjoint  operator  matrix  equation  (49)  have  the 
following  relationships; 


<t*u1,j-'5j>R  *  -  <*'“(•  tij.j’R 


+  tlj>S1  *  <l*“f  Ifjnj>S4 

*  <t*,u("J»'-rtj>s11  +  fr-,j"j)'>s4.  <56 ) 


<l*V  1>R 


*  <t*wjni,77>s2  ♦  <t*w1,jrn,>S3 

I  I 

+  <t*(w.n. )  ,TT>n  +  <t*w. ,  ( 7r  n  . )  V 

l  l  52i  i  i  i3i 


(57) 


In  writing  (56)  and  (57)  we  assume  that  <  ,  >R  can  be  evaluated  as  the 
sum  of  quantities  evaluated  over  subregions  of  R  such  that  all  the  sur¬ 
faces  S^. ,  S2-,-,  S3l-,  S^-  are  contained  in  the  union  of  the  boundaries  of 
these  subregions. 


4.3.2  A  General  Variational  Principle 


For  the  operator  equations  (49),  we  define  the  governing  function, 
following  (45)  as; 


1 


0(u)  =  <Pu,,u,>R  ♦  2< #>2wi * U1  >R  -  ^tjjj-VR  *  «-  *  “i>, 


-  ^ i ,wi >R  +  <t*'*M'ff>R  -  2<t*  f’n>R  * 


■  2<t’eij-rij>R  +  <t*(E,jk1  ♦  a  M  8(j  -6kl)  ekl,  e(j>R 


♦  2<t*a  M  Sjjejj.f  >R  +  <t*M£,f>R 


-  2<Wr  -  2<Wf 


-  <  ^■j-t*(ui-2ui)nj>S1  -  <».t*(wj-2wf)nf>S; 


*  <V  t«(Jr-2Jr>v$3  .  <Uf,  t*(  TunJ-2T,)>S4 


<riJ,t*((ujnjr-2(giljnj)>si.  -  <W.t*«(w1nf  )-292>>s 


2i 


+  <w. ,  t*((  7rn.  )’-2g  n.)><.  +  <u-,  t*((  X<  ini  >  ‘-29  "<)><; 

1  1  3  1  53i  1  J  4  1  54i 


(58) 


The  Gateaux  differential  of  this  function  along  v  =  [u^ .w^ ,  7T,^j  . ,e.j  . ,^j 
is; 


\aM  -  <0,.,  pUj  *  p2w,  -  -  2Ff>R 


*  <uf  *  p2" i  -  ‘‘rij./R 


I 


*  ‘V'V'i  *  "  l*  Tf  ■  2Vr 

P-  1 


*  <»,,»«.  *  ( -  +  1* - )w.  -  t*  7T  (>n 


v.1 


*  <7r ■  ‘**1,1  -  *  <V-  -  ‘*f>R 

+  <>Eij>  t,uij  -  ‘‘‘u'r  ♦  <rij-  -  t*Vn 

*  <*1.1-  -‘*^1  *  W£.m*  «  HSfAilekl  +  aM8<if 


i  j  ^  R 


+  <6ij‘  ’**ii  +  «  M5ij5kl^kl  +  ««*iJ?>R 

+  <?,  -t*7T  +  t*  aM  5klek1  +  t*  m£>r 

+  <£  ,  -t*7 r+  t*  a  M  6k]ekl  +  t*  m£ >R 

-  'V  ‘‘'“("j  -  ^V’s,  -  <r,j-  t*Il"j>S1 

-  <*,  t*(V,  -  Zw^O^  -  <JT,  t*w,nj>S2 

+  <W.#  t*(7rn.  -  27Tn.j )>^  +  <Wj,  t*7rn.>s 

3  3 

<nf  -  2V>s4  +  <uf  t*rijnj>s4 

-  <rfJ,  -  <tu,  t*(u)nJ)'>3ij 

-  <7f,  t*(  (w.n. ) 1  -2g  ><•  -  <7 T,  t*(w.n.)'>c 

i  1  2  b2i  11  b2i 

+  <w-,  t*((7rn. ) '  -2g  n.)>r  +  <w, ,  t* ( 7r  n  - )  * > «- 

1  1  3  1  b3i  1  1  b3i 

+  <15,- ,  t*(('Ci1ni),-2(g  )n.)>,  +  <u.,  t*(t\Mn,)  '>.  (59) 

1  'j  J  4  '  '>4i  1  J  i4j 


Using  equations  (56)  and  (57),  the  Gateaux  differential  can  be  rewritten 


\n,u)  =  2<uc  (,ui  *  Vi  -  ‘"Cfjj  -  Ff 


I 


R 


i,  <*,4*^' 


i r 


ff 

i’ll** 

$ 
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_  ^2  1 

+  2<w. ,  Pu,.  +  (— =  +  l*~)w.  -  t*  7T,  -  G*>„ 

II  ^  k  ^  >  i  i  k 

+  2<if-  <**1.1  -  >R  *  2<*1J-  ‘’“i.j  •  **ei j’r 

*  2<ftj-  -‘*^j  *  0  MSiAi)eki  *  oMM  >r 

+  2<£ ,  -t*7 r+  t*  a  M  6k1ekl  +  t*  m£  >R 

-  2<ttj,  t*(Ujnj  -  u,nj)>Si 

-  2<W  ,  t*(w.n.  -  w.n. )>c 

li  ii  ^2 

+  2 <w . ,  t*(7Tn,  -7Tn.  )>c 

*2<V  flTijnj  -t,)>s4 

-  2<rtJ,  t*((u1nJ)'-(gi>1nJ)>sii 

-  2<7 r,  t*( (w.n. ) '-g  >«. 

1  1  2  b2i 

+  2<w. ,  t*((7rn.)'-g  n.)><- 
1  1  3  1  ^3i 

+  2<u.,  t*((  r-  .nJ’-fg  )n,)>s  U 

'  ij  J  4  1  ->41- 


The  Gateaux  differential  vanishes  if  and  only  if  all  the  field  equations 
along  with  the  boundary  conditions  \ 54 )  and  the  jump  conditions  (55)  are 
satisfied  because  of  linearity  and  nondegeneracy  of  bilinear  mapping 
(46).  Hence,  vanishing  ofAp(v)  for  all  v£W  implies  satisfaction  of 
(49),  (54)  and  (55). 
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4.3.3  Extended  Variational  Principles 

Equations  (56)  and  (57)  relate  pairs  of  operators  in  the  operator  ma¬ 
trix  (49).  These  relations  may  be  used  to  eliminate  either  of  the  ele¬ 
ments  in  each  pair  of  the  function  fl(u)  in  (58).  Eight  alternative  forms 
can  be  obtained  by  using  either  or  both  relations.  Elimination  of  an  op¬ 
erator  A.,  from  the  function  implies  that  state  of  variable  u.  needs  not 

I  J  J 

be  in  the  domain  M . •  of  A-..  Where  A.  •  are  differential  operators,  this 

'  J  *  J  '  J 

results  in  relaxing  the  requirement  of  smoothness  in  u-  thereby  extend- 

•  J 

ing  the  space  of  admissible  states.  In  the  context  of  finite  element 
method,  it  is  clear  that  the  extension  of  the  admissible  space  provides 
greater  freedom  in  selection  of  approximation  function.  In  the  follow¬ 
ing,  the  possible  extensions  are  explicitly  stated. 


Using  (56)  to  eliminate  ,  from  (58), 

*  J  >  J 


ni  ’  <  "“i-Vr  *  2<p*VVr  +  <(-  *  I'-Wr 

-  <t*n>VR  +  <t*"l,i-,r>R  -  2<t*f -^R  *  2<t*ui,j-'Clj>R 

-  2<t*efj-riJ>R  +  <t*(Eijk1  +  a  M  6ij  Skl  >  ek1-  e1  j>R 
+  2<t*a  M  jei j » f  >R  +  <t*M£,£  >R 

-  2<u.,Fi>R  -  2<w.,Gi>R 


-  2 


<rij’t*(ur°1)nj>S1  -  <”’,t*(wr2fi,.)ni>S; 


+  <V  t*(!r-2ir)nf>S3  -  2<Ul,  t*Tj>s^ 

-  Z<H. . , t*( (u -  n  • ) ' - ( g  )  - n - ) > o  -  <7r,t*((w-n, )  '-2g  )>, 

U  i  J  i  1  J  bii  11  2  b2i 


*  V  *  '•-  *  V  ' 


(61) 


+  <w.,  t*((7rn.  )'-2g  n.)><.  -  2<u.,  t*g  n.><. 

1  1  3  1  b3i  1  4  1  b4i 

In  (61)  the  stress  components  need  not  be  differentiable.  Using  (56)  to 
eliminate,  alternatively,  u,  .,  (58)  gives 

*  9  J 


Oo  1 

n2  =  ^Wr  +  -  2<t*'rij,j-“l>R  +  <(-  +  1*->Wr 

-  <l*  "i-Vr  +  <t*“lJ-,r>R  ‘ 

-  2<t*efj-'tij>R  *  <l*(E1jkl  +  “  H  51j  6kl>  *k1>  Vr 

*  2<t*«M  f>R  ♦  <»*"*.  f>R 

-  2<Wr  -  2<1VVr 

*  2<'fiJ.t*Vj>s1  - 

+  <V  t*(7r-2fr)„.>$3  +  2<u.,  t*<t(jVV>s4 

+  2<t.-,-,t*{g  )  -n  •><-  -  <7r,t*((w.n.)‘-2g  )>«. 

1J  l  1  J  11  2  b2i 

+  <w.,  t*((7rn- )'-2g  n.)><-  +  2<u- ,  t*(  ( •  -n  - )  ‘  ~g  n-)><.  (62) 

i  i  3  1  ^3i  1  U  J  4  i  ^4i 

in  which  the  displacement  need  not  be  differentiable.  Elimination  of 

w.  •  by  (57)  from  (58)  gives 
1,1 


n3  *  ^i-Vr  *  2<»»f.«i>R  -  ‘‘"t'ijj-VR  *  «~f  +  r->W 

-  2<t*  ”>VR  -  2<t*f  •’r>R  + 

-  2<t*etj-rij>R  *  ^ijkl  *  0  M  Sfj  5kl>  ekl’  ei j>R 


.  V  -  X  -  j.  -V  W  J 


+  2<t*o.MSjjejj,f>R.<t*Hf,f>R 

-  2<u,,F1>r  -  2<w1,G(>r 

'  <tfj-t*lur2°1)nj>S1  +  2<’r"t*>'?nf>S2 
+  2<V  t*(  IT- jr)n(>s^  ♦  <«,,  t*(T,Jnj-2t1)>s4 

-  <rii.t*({u.n.),-2(g  ).-n.)>,  +  2<7r,t*g  ><- 

■  J  i  j  i  '  J  ^ 1 i  2  ^2i 

+  2<w. ,  t*((7rn.)’-g  n.)>-  +  <u_-,  t*( ( 't'.- ,-n  - )  * -2g  n.)>r  (63) 

l  1  3  1  b3i  1  J  J  4  1  54i 

Here,  w.  need  not  be  differentiable.  In  the  same  way  tt  •  can  be  dropped 
■  >  • 

out  by  using  (57),  yielding 

P?  1 

fl4  ■  ^WR  *  2<^2W1  *ui>R  -  +  <(~  +  1*-)W,.»,>R 

+  2<t*wi,i-,r>R  -  2<t*f.>F>R  *  <t*u,j.t,j>R 

2 

-  2<t*e(j,t,j>R  *  <t*(EMk1  *  a  N  «„  5kl)  e,,.  ei(>0 


ij  kl '  kl  ’  1  j  R 


*  +  >R 

■  2<Wr  -  2<Wr 

-  <r)J,t*(uj-2tilnj>Si  -  2<W,t*(«1-S1)n1>S2 

-  2<V  ‘**Vs3  *  <V  t*<tu"j-2fi,>s4 

-  <T1-1-,t*((u.n.)’-2(g  )  .n .  )>(.  -  2<  7T,t*(  (w.n, )  '-g  )><- 

ij  i  j  1  *  J  ^l*f  11  2  **2i 

-2  <w . ,  t*g  n  .>s  +  <u • ,  t*( (  X*  ) *-2g  ”•)>, 

i  3  1  b3i  1  J  4  1  b4i 


ffr  t?  ikrTin .  rfr  ^  ^  "  E  -  I 


»»?  •**  *»W  *IF 


In  (64)  the  fluid  pressure  is  not  required  to  be  differentiable.  As  can 
be  seen  from  (61)  to  (64),  both  the  differential  operators  in  an  adjoint 
pair  cannot  be  removed  at  the  same  time.  Use  of  (56)  and  (57),  however, 
eliminates  the  differentiability  of  two  field  variables  from  (58).  Elim¬ 
inating  T,-.-  .•  and  tt  .•  from  (58),  we  have 


p2  1 


n5  *  <p  Wr  +  2<|,,Wr  +  <(~  +  i*T)wi’wi 


.->r 


♦  2<t*wf>i,ir>R  -  2<t*^,7T>R  ♦  2<t*UlJ,Tfj>R 

-2<t*ejJ.,  r(J>R  +  <t*^i  jiti+  5,  j  Skl '  ekl'  ei j'r 

*  2<t*aH  «(JefJ,f>R  *  >R 

-  -  2<w1-g1>R 

-  2<T1j,t*(u1-01)nJ>si  -  2<7T,t*(wr«,)n,>S2 

-  2<w(,  t*#nf>Sj  -  2<uf ,  t*T()>Sii 

-  2  <  T . . ,  t*  ( ( u  •  n . ) '  -  ( g  ).n.)><.  -  2<7T  ,t*(  (w-n, )  '-g  )><- 

1  j  'j  1  ^  i  i  2^2 i 

-  2<w . ,  t*g  n  ■><•  -  2<u . ,  t*g  n.)><. 

1  3  1  b3i  1  4  1  b4i 


(65) 


In  (65),  the  total  stress  field  and  the  fluid  pressure  need  not  be  diff 
erentiable.  Elimination  of  u.  .  and  IT  ■  from  (58)  gives 


n6  *  <*,ui*u1>R  *  2<<’«f*u1>R  -  <t*'Cjj,j.u,>R  *  <(-  ♦  !*-»», ."i>R 


*  -  2<t*M7->R 


-  2<t*efj-  "Vr  *  ^ijkl  *  «  M  Sij  Skl>  ekl-  ei 


*2<l*0,lVii'f,«  *  <t*M  f-f>R 


-  2<VVr  -  2<«1-s1>R 

*  2<T1J,t*u>j>s1  ■  ZOT.tMw,-#,)!,,^ 

-  2<wf>  t*jrn,>$3  *  2<u,,  t*(r,j"j-T,)>S4 

*  2<riJ,t*(gi)jrj>Sij  -  2<7r,t*((w1.n1)'-g2)>S2f 

-  2<w,,  t*g3n,>S3j  *  2<u,, 


4  1  4i 


where  ui  and  ir  need  not  be  differentiable.  Eliminating  X 
from  (58), 

P2  1 

n7  =  +  2<Pzwi-ui>R  *  <(~  +  !*-)»( .*»|>r 

-  2<t*  -  2<‘*f."r>R  *  2<t*“f,j.'C1j>R 

'  2<t*eij'Tij>R  +  “M  5ij  5kl>  ekl-  eU> 

♦  2<t*a  M  >R  *  <t*M  f.f  >R 

-  2<Wr  -  2<Wr 

-  2<'Cij"'*,“i'u'i>nj>S1  *  2<Jr-tSni>S2 

*  2<",,  -  2<u,,  t*f j>$4 

-  2<T-i,t*((u.n.),-(g  )  .n  .  )><-  +  2<7r,t*g  ><- 

ij  i  J  i  1  J  hi  2  b2i 


W 


ft  \ 


+  2<w. ,  t*((7rn, )‘-g  n.)>r  -  2<u.,  t*g  n,)>«.  (67) 
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which  does  not  require  the  total  stress  and  the  relative  displacement  of 
fluid  to  be  differentiable.  Using  (56)  and  (57)  to  eliminate  u.  .  and 

'  f  J 

w.  (58)  is 

'  J  * 

1 

n8  ■  <puj »ui>r  *  ^VVr  -  +  <(-  * 

-  <t*irt1.wf>R  -  2<t*f  ,7t>r 

-  2<t*ei j-'rij>R  *  <,;*(E1jkl  *  «  M  6ij  5kl>  'kl-  *fJ>R 

’  2<u1  *E1>R  -  2<1Wr 

*  2<  1j,t*u?"j>S1  * 

*  2<V  )n(>$3  -  2<V 

+  2<Tii»t*(g  ).-n  •><•  +  2< 7r,t*g  >, 

ij  j  i  J  i>li  2  52i 

+  2<w. ,  t*(  (  7rn. ) 1  -g  n.><-  +  2<u,,  t*(  (  Zi  -n  .)  ‘-g  n.)><-  (68) 

l  1  3  1  *31  1  iJ  J  4  i  i4i 

Here,  the  solid  displacement  and  the  relative  displacement  of  fluid  need 
not  be  differentiable. 
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4.3.3  Specializations 

If  the  admissible  state  is  constrained  to  satisfy  some  field  equa¬ 
tions  and/or  the  boundary  conditions,  certain  specialized  forms  of  the 
variational  principle  are  realized.  This  procedure  is  used  to  reduce  the 
number  of  free  variables  in  the  governing  function.  Also,  certain  as¬ 
sumptions  in  the  spatial  or  temporal  variation  of  some  of  the  variables 
lead  to  approximate  theories.  In  the  context  of  direct  methods  of  ap¬ 
proximation  the  constraints  assumed  in  the  specialization  must  be  satis-- 
fied  by  admissible  states.  If  it  is  difficult  to  satisfy  the  con¬ 
straints,  such  specialization  of  the  variational  formulation  will  not  be 
useful  in  practice.  Some  specializations  of  the  extended  variational 
stated  in  the  previous  section  are  presented  below. 


For  the  functional  in  Equation  (65),  in  which  the  soil  stress  and 
the  fluid  pressure  need  not  be  differentiable,  specialization  to  satisfy 
( 49 ) 3  and  ( 49 )^ ,  i.e.  satisfying  identically  the  kinematic  relationships 
gi  ves 


0,  =  <Pu(,u,>R  +  2<e2  VVr  +  <(-  *  l*-)wfVR 

+  o  M  5ij  6kl>  ekl  ’  e1j>R 


+  2<t*a  M  5. ,e, ■ 

'  J  '  J 


£>r  +  <t*M  >R 


"  2<Ui »Fi >R  -  2<Wi»Gi>R 

-  2<t,J.t*(ur8,.)nj>Si  -  2<ir,t*(»r«,.)nf>S2 

-  2<w. ,  t*£n.>c  -  2<u  • ,  t*f.)><; 

1  i  ^3  1 


-  2<  r.  .,t*((u-n.)  -(g  )  -  n  •  )>«-  -  2<  n  ,t*(  (w.n . )  *-g  )><. 

ij  1  J  i  >  J  ^ 1 i  11  2  ^2i 


2<w.,  t*g  n  •>«.  -  2<u . ,  t*g  n.><- 

1  3  1  b3i  1  4  1  b4i 


If  the  field  variables  over  the  domain  are  continuous,  the  jump  discon¬ 
tinuity  terms  drop  out  giving  the  specialization; 


_  p?  1 

%  ■  <pWr  *  2<p2  Wr  *  «“  +  >*->Wr 

+  <,;*(Eijkl  +  aM  4kl>  «kl-  e1j>R 
*  2<t*aM  5ijeij>f>R  + 


-  2<Uj-F(>R  •  2<Wj.  Gi>R 

-  2<rij-t*,uf-0i,nj>s1  -  2<*.t*(wr«()iySi 

-  2<w,,  -  2<Uf,  t*T()>$4 


Further  specialization  of  (70)  to  the  case  where  displacement  boundary 
conditions  are  identically  satisfied  yields  the  function  governing  the 
two  field  formulation  proposed  by  Ghaboussi  and  Wilson  [12]  except  that 
in  the  present  formulation  the  boundary  terms  are  consistent. 

Alternatively,  specializing  Equation  (67)  to  satisfy  the  (49)4, 


nn  ■  <pWr  *  2<p2  Wr  +  «-  *  i*— >R 

'  2<t*7r,Cwi>R  *  2<t*f.”->R 


*  <t;*,Eijkl*  aH  6ij  5kl>  ekV  e1j>R 


'»SN 


*  2<t*aMVlj'f>R  + 

-  2<u,,F(>r  -  2<w,,G1>r 

-  2<-rfj,t*(urtij)nj>Si  ♦  2<7r,t*vvyS2 
+  2<w1  ’  fUr-frln,^  -  2<Ui>  t*f,>S4 

-  2<T1-i,t*((u.n.),-(g  )  -  n  .)><-  +  2<7T,t*g  ><• 

U  i  j  l  1  J  bii  2  b2i 

+  2<w. ,  t*( (  7T n  - ) '  -g  n.  )><•  -  2<u.,  t*g  n.)k  (71) 
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Furthermore,  assuming  that  the  internal  discontinuities  and  the  boundary 
conditions  on  Sj  and  are  identically  satisfied  and  eliminating  £  by 
using  (49)g,  gives 


1 


ni2  =  <  PU1-UT>R  +  2<P2  Wr  +  <(“  +  l*-)wf,wf>R 


-  2<t*  VVr  +  <l*2i jklekl ’  «1j>R  +  2<t*“51j  e1j-7r>F 


-  <t*7T/M,!r>R  -  2<uj,F1>R  -  2<wj,G1>, 


*  2<*'-t*wini>S2  '  2<V  t’Vs, 


(72) 


which  is  three  field  formulation,  has  important  meaning  in  the 
context  of  finite  element  analysis,  since  thismixed  formulation  of  u-w- 
V  can  produce  the  continuity  of  the  pore  pressure  which  is  the  physical¬ 
ly  important  quantity  in  the  analysis  of  dynamic  response  of  the  fluid- 
saturated  solid.  Similar  three  field  formulation  can  be  obtained  by 
specializing  n,  to  satisfy  (49k,  (49).  and  (49k. 


Clearly,  a  large 
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number  of  other  specializations  are  possible  even  if  they  are  not  listed 
here  and  left  to  interested  readers. 


4.4  COMPLEMENTARY  FORMULATION 

4.4.1  Complementary  Form  of  Field  Equations 

An  alternative  procedure  to  set  up  the  variational  principles  govern¬ 
ing  the  problem  is  to  write  the  operator  equations  in  complementary  form 
instead  of  the  direct  form  of  equation  (49).  In  writing  the  equations, 
it  is  assumed  that  kinematics  of  the  solid  and  the  fluid,  i.e.  (30)  and 
(31)  are  satisfied.  The  complementary  form  of  field  equations  is 
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(73) 


where 


P  =  -t*(l/M  ♦  or  Cijkl5kl5ij) 

L  =  —  t*(5.m  +  5.  2L) 

2  1  aj  Jm  ai 


(74) 

(75) 


The  governing  function  for  the  set  of  equations  (73)  along  with  the 
boundary  conditions  and  internal  jump  discontinuity  conditions,  follow¬ 
ing  Sandhu  [11],  is 


J  » 


<PWr  +  2<  P2  Wr  -  “i 


,,P2 


*  <(-  *  1*-)W(,W(>R  -  <t*  ir,,w(>R  *  <t*wM,77>R 


■  <t,(  If  *  ^  C1  jklSk1S1  j )7r  ■  ^  'r 


*  2<t*“  Cf  jkl5kl  1j-7r>R  +  <t*u(J-'rfj>F 


‘  <t*Cijkl  kl’l!  "  2<ui,Fi>R  "  2<wi,Gi>R 

'  <rij't*,ui'2<!ilnj>S1  -  <’1'.t*(w1-2ft1)n1>S2 

♦  <«,,  rlir-ur  In,^  ♦  <u,, 

-  <Tij't*,,ufnj)'-2(91)fnj)>Su  -  <7T.t*((wfnf  )--292>>s 


2i 

+  <w.,  t*(  (  7rn  - )  *-2g  n.)><.  +  <u.,  t*(  ( T-  ,-n , ) 1  -2g  n,  )>«-  (76) 
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As  in  direct  formulation,  it  can  be  shown  that  the  Gateaux  differential 
of  (76)  vanishes  if  and  only  if  the  field  equations  (72),  the  boundary 
conditons  (54)  and  the  jump  condtions  (55). 


4.4.2  Extended  Complementary  Variational  Principles 

Following  the  principles  and  methodology  presented  in  Section  4.3,  it 
is  possible  to  develop  extended  variational  principles  for  the  comple¬ 
mentary  form  of  the  field  equations  (73)  as  well.  Relations  (56)  and 
(57)  can  be  used  to  eliminate  some  of  the  operators  from  (76).  As  a  re¬ 
sult,  following  extensions  of  (76)  are  possible. 

Use  of  (56)  an  (57)  to  el  iminate  T- ^  ,•  and  77”, -from  (76)  qives 


J;  =  <Puj,ui>R  ♦  2<  P2  ♦  <(— * 

*  2<t*wM,7T>  -  a2  Cljk1Skl«1J)W.7T>R 

*  2<t*a  Cijkl5kl  fj-7r>R  +  2<t*ui,j.'Tfj>R 

"  <t*Cijkl  ’  2<ui  »Fi>R  "  2<wi  >Gi  >R 

-  2<ri.j,t*(ui-01.)nj>Si  -  2<7r,t*(wrwi)n1>S2 

-  2<w.,  t*7Tni>s^  -  2<Ui,  t*VS4 

-  2<T.  .,t*((u.n  •)  *-(g  ),ni)><-  -  2<  7T,t*(  (w,n, )  ‘-g  )>«. 

ij  ij  i  i  j  *■>  1 1*  2  2i 

-  2<w. ,  t*g  n,>c  -  2<u.,  t*g  n->c  (77) 

1  3  1  oi  1  4  1  b4i 

Here,  the  stress  field  and  the  pore  pressure  need  not  be  differentiable. 

Similarly,  elimination  of  X..  ■  .  and  w.  .  from  (76)  leads  to 

I  J  >  J  I  J  • 


J2  =  <  PO, ,uf>R  *  2< P2  »,,U(>R  *  «-♦  r- )«,.»,>„ 
-  2<  t*  7r  j  ,w,  >  -  £  CiJkl5kl5ij),r',r>R 

*  2<t*aCijkl5klT1j>’T'>R  +  2<t*“fJ-^j>R 


'  <t*Cijkl  ^"ij’^kl  >R  '  2<ui'Fi>R  '  2<"i'Gi>R 

-  2<TfJ,t-(urC,-)nj»Si  *  2<7r,t*;jni>S2 

+  2<w . ,  t*(  7T-1T  ) n .  > „  -  2 < u  ■ ,  t*T  . >c 

1  1  I  ^ 

-  2<T.1-,t*((u.n.),-2(g  ).n.)>«.  +  2<7T,t*a 

ij  i  J  1  1  J  b 1 i  2  b2i 


■  <t:*(~+  0  Cijkl5kl5ij)7r>7r>R 

+  2<t*aCijkAl  V^R 

<t*C i jkl  j ’ ^*k 1 >R  "  2<ui*Fi>R  '  2<wi’Gi>R 


+  2<rij*t*u?nj>s1  +  2<7r*tSni>s2 
+  2<V  t*(7T-^f)n.>S3  +  2<uf,  t*(TfjnrTf)>s4 

+  2<t.-,,t*(g  ).-n.>,  +  2<7T  ,t*g  >, 

ij  i  1  J  2  *2f 

+  2<w. ,  t*((7rn.)'-g  +  2<u-,  t*((T.- )  *-g  n.)><-  (80) 

1  1  3  53i  1  iJ  J  4  i  ^>4l- 

Derivatives  of  ik  and  are  not  involved  in  (80). 

4.4.3  Specializations  of  Complementary  Extended  Variational  Principles 
As  in  direct  formulation,  some  specializations  on  the  complementary 
extended  variational  principles  are  possible  by  requiring  that  certain 
field  equations  and/or  boudary  conditions  and  jump  conditions  be  identi¬ 
cally  satisfied.  Here,  we  present  some  secializations. 

If  we  assume  that  (73)4  is  identically  satisfied,  Jj  results  in 

p2  1 

J5  -  <  pu,,vR  *  2<p2  Wr  *  <(-  *  ‘'-'Wr 

*  -  <t*(-jU  a2  Ctjk1*k1«1J>W.W»R 

+  <t*C1jklflj.t'kl>R  -  2<«,.F|>r  -  2<"i.Vr 
-  2<TfJ,t*(Uj-ttf)VSi  -  2<7r,t*{wrw,.)n(>S2 
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-  2<w,.,  t*#nf>S3  -  2<u,, 

-  2<Tii,t*((u,n-)'-(g1  )<ni)>c  -  2<ff,t*((w,n.)'-g  )><- 

ij  ij  1  1  ^  ^li  ’  1  2  °2i 


2<wi 


t*g  n,>s 
3  1  b3i 


-  2<u • ,  t*g  n,>s 
1  4  1  b4i 


(81) 


Specializing  Jj  to  satisfy  ( 73 ) 3 , 


J6  -  <  Puf,u,>R  *  2<p2  w,,Vr  *  «-2  ♦  V*-i)wj,w1>R 

*  <t*<-U  a2  +  2<t*unJ,  r(J>R 

-  <t*CiJkl‘fl4*rkl>R  -  -  ^VVr 

-  2<Tjj,t*(ur0i)nJ>Si  -  2<Jr,t*(wi-S1)iyS2 

-  2<wf,  finy^  -  2<u,,  t*ff>$4 

-  2<rjJ,t*((ujnJ)'-(g1))nJ.)>si.  *  2< "\t*< («(",) 

-  2<w,,  t*g3P,>s3.  -  2<uf.  t*34"i>s41 


(82) 


It  is  unlikely  from  the  above  two  specializations  that  further  speciali¬ 
zation  on  Jg  or  Jg  through  requiring  the  other  field  equations  to  be  sa¬ 
tisfied  would  make  them  simpler.  Since,  moreover,  specializations  of  the 
extended  functionals,  to  would  show  little  differences,  we  do  not 
present  every  possible  specializations. 


SECTION  V 


DISCUSSION 

A  systematic  development  of  variational  principles  for  linear  elasto- 
dynamics  of  fluid-saturated  solids  have  been  presented.  Nonhomogeneous 
boundary  conditions  and  internal  jump  discontinuities  were,  explicitly 
incorporated  in  general  variational  principles.  Allowance  of  jump  dis¬ 
continuity  terms  in  variational  formulation  is  meaningful  in  the  context 
of  direct  approximation  in  finite  element  spaces,  since  the  space  of  ap- 
proximants  may  not  be  sufficiently  smooth.  Based  on  the  direct  and  the 
complementary  formulations,  extensions  of  the  variational  principles 
through  elimination  of  certain  field  operators  and  specializations  by 
restricting  some  of  the  field  equations  and  boundary  conditions  to  be 
identically  satisfied  have  been  proposed.  These  formulations  should 
provide  a  basis  for  development  of  approximate  solution  procedurse  and 
also  approximation  theories  governing  the  problem.  Fig.l  diagramatical- 
ly  depicts  the  possible  extensions  of  the  general  variational  principle 
based  on  the  direct  formulation.  Fig. 2  shows  the  same  for  the  comple¬ 
mentary  formulation.  In  either  case  the  specializations  listed  in  this 
reportare  shown.  Evidently,  other  extended  forms  could  be  used  as  start¬ 
ing  points  for  specialization.  In  the  interests  of  brevity  we  desist 
from  attempting  to  catalog  all  the  possibilities. 


Satisfy  u,w  B.C. 
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